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A Gaussian variational approximation is often used to study interfaces in random media. By 
considering the 1+1 dimensional directed polymer in a random medium, it is shown here that the 
variational Ansatz typically leads to a negative variance of the free energy. The situation improves by 
taking into account more and more steps of replica symmetry breaking. For infinite order breaking 
£^ ' the variance is zero (i.e. subextensive). 

This situation is reminiscent of the negative entropies in mean field spin glass models, which were 
^ ■ also eliminated by considering infinite order replica symmetry breaking. 
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The subject of interfaces in random media has many application, such as directed polymers in a random medium, 
domain walls is dirty magnets, and interfaces in random sponges. Various aspects of this field and methods of approach 
have been reviewed, see Q, ||. 

As a tool to get a grip on these complicated objects, a Gaussian variational Ansatz has been proposed by Mezard and 
Parisi Q to calculate the free energy. It was found that replica symmetry is broken spontaneously. This variational 
replica-approach reproduces the Flory values for critical exponents, that can be derived by Flory-Imry-Ma type of 
■ estimates. 

In the approach an infinity (fc = oo) of replica symmetry breakings occur. Carlucci, de Dominicis and Tcmesvari 
have shown that the instability of the replica approach at finite k disappear in the limit k — > oo. [^) 

The problem of a directed polymer on a random substrate in d = 1 + 1 has received quite some attention. It was 
pointed out by Kardar S that the replicated free energy can be related to the ground state energy of n interacting 
bosons in one dimension. This problem can be solved exactly by means of the Bethe Ansatz. The depinning of such 
a polymer or interface from a wall is also solvable exactly, when the effect of the wall potential is described by a 
derivative condition of the wavefunction at the wall. Various interesting results have been derived from the exact 
' groundstate wavefunction near the depinning transition. 

Here we reconsider the d = 1 + 1 problem of a directed polymer of length L in a random potential. The purpose 
is not to derive the best variational approximation to the exact solution. Having in mind more complicated models, 
we wish to test a particular aspect of the variational approach of Mezard and Parisi, that we fully work out for 
this situation. Due to the one-dimensional character of the problem, we can get the results by studying a quantum 
k> , Hamiltonian. 

In particular we point out that order n L terms of the n-fold replicated free energy F n tend to have the wrong 
sign, so that they tend to predict a negative variance of the physical free energy. This unphysical prediction occurs in 
case of no replica symmetry breaking and becomes less severe for a finite number of breakings; it seems related to the 
instability of some fluctuation modes. The whole n 2 L term disappears if infinite order replica symmetry breaking is 
taken into account. The negative variance of the free energy thus plays a role similar to the negative entropy of the 
SK- model: both vanish only when infinite replica symmetry breaking is taken into account. 

The implication of vanishing of the n 2 L term is that the fluctuations in the free energy are smaller than L 1 / 2 . 
We shall then show that the replicated free energy has leading correction ~ n 5 L, which means that the free energy 
fluctuations are of predicted to be of order L 1 / 5 . This is in agreement with the Flory prediction worked out by Mezard 
and Parisi. Note that the exact order of magnitude is known to be L 1 / 3 . 
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In section 1 we introduce the model and the variational Ansatz. The simplest cases are considered. In section 2 
the problem is solved for infinite order replica symmertry breaking. We close with a summary in section 3. 



I. THE MODEL AND THE VARIATIONAL APPROACH 

We consider an interface z(x) with Hamiltonian 

f3H = £ dx^^j +V(x,z(x))} (1) 

(3 — l/T is the inverse temperature, 7(T) is the interface stiffness, and there is 5-correlated Gaussian disorder 

V{x, z)V{x', z 1 ) = <7 2 5{x - x')5{z - z') (2) 
where the overline denotes the average over the quenched disorder. We shall study the partition sum 

Z = Jdz e-W = e-? F (3) 

(For a formulation on the lattice, see ||). The replicated partition sum 

oc 

Z n = Z* — » exp(-/3F„) = exp(-n/3F) = exp£n fc (-/3F) c fe um (4) 

k=l 

generates the cumulants of (3F as the expansion coefficients of j3F n for small n. For having a non-negative variance 
of F the n 2 term of F n should be non-positive. 

Averaging over disorder one obtains the effective n-particle Hamiltonian 

n n = / dx ii E(^) 2 " ° 2 £ - z fj {x))} (5) 

By analog with the Feynman-Kac path integral formulation of quantum mechanics, the replicated free energy can 
be obtained as 

(3F n = LE n (6) 
where E n is the ground state energy of the quantum Hamiltonian 

I n d 2 

' a=l a a<p 

The expansion 

E n = nE<V + n 2 E^ + n 3 E^ + n*E^ + n 5 E^ + ■■■ (8) 

determines the cumulants 



(0F)L m = {-if-'klLE^ (9) 
The present problem has an exact ground state wavefunction 



ipo(z) = exp-K £ \z a - zp\ (10) 

a<(3 



where k — a 2 j. The groundstate energy is 



^ = ^-f (ID 
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To our knowledge, the variational approach to this problem was first formulated by Honeycutt and Thirumalai [Q. 
They considered a Hartree variational interface Hamitonian, that we generalize as 



r L " /h \ 2 " 

• /U o=l V 7 a,/3=l 



(12) 

a,/3=l 

The original approach had q a p — 8 a pq. The generalized approach is equivalent with the one of Mezard and Parisi, 
though it is simpler. Indeed, the mapping onto a quantum problem, which can only be done in d = 1 + 1, prevents 
the introduction of a Fourier wavenumber. 

The variational Hamiltonian is equivalent with considering the Hartree variational quantum Hamitonian 



2 7 ^ dzl 2 ^ 

which has the groundstate wavcfunction 



= -^E^-?E »«(« 2 )«.w» (13) 



The ground state energy now is approximated as 

S„=^ + (^1(^-^)1^) (15) 



Rescaling q — ► 2er 4 7 2 /7r we have 



with 
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= (16) 

Z7T 



Bn = T q aa - Y (17) 
The exact value would correspond to E n = 7r(n — n 3 )/3. We shall now consider various choices for the shape of q a p. 



A. Diagonal q 

The case q a p = qd,S a ,i3 was studied by Honeycutt and Thirumalai. Optimizing E n yields q d = (n — l) 2 , so that 

E n = -n{n-lf (18) 

which differs in sign from the Honeycutt-Thirumalai result. Since our n 2 term is positive, we conclude that this result 
leads to a negative variance of F, and is not acceptable. 



B. Translational invariance and off-diagonal q 

1. Replica symmetry 

If there is no RSB then we can set q a p = (q d — qa)5a,p + <Zo- 

E„ = nq d - 2n(n - l)^q d - q (19) 

The structure of this expression shows that we only get a well posed problem if we relate go and qd in some way. 
Above we considered qo — but the result was unsatisfactory. As we have averaged over disorder, there is translational 
invariance. One would expect that the physically relevant variational Ansatz reflects the translational invariance of 
the replicated problem. Hereto we have to impose 
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^ la/3 = (20) 


This condition is equivalent to staying in the "replicon" or "ergodon" subspace. 

implying here that qa + (n — l)qo = 0. This allows to eliminate qd, after which we get 

E n — —n(n - l)q - 2n(n - l)y/-nq (21) 

The saddle point is go = —n, yielding 

E n = -n 2 (n-l) (22) 

This result has E = and variance —1, in both respects better than previous case. But also this result leads to a 
negative variance, and is not acceptable. 



2. One step replica symmetry breaking 

Now we choose 

q a p = (qd - qi)S a fj + (qi - qo)?^ + qo£%p (23) 

where £( Xl *> is the Parisi 1RSB matrix having n/x\ diagonal x\ x x\ blocks with elements equal to unity. (Thus £W 
is the identity matrix, while is a matrix with all elements equal to unity.) Due to translational invariance we 
may eliminate qd and get 

qd = -{xi-l)qi-(n-x 1 )q (24) 

The inverse reads 

-1 _ <W ffl - gQ £(xi) 90 g(n) 

a(i qd-qi {qd - Qi){Qd - 0- - xi)qi - Xiq ) afi (qd - (1 - X\)q\ ~ Xiq )(q d - (1 - Xi)qi - x x q + nq ) afi 

(25) 

The last factor diverges in case of translational invariance. Fortunately, it generally cancels from our expression for 
E n . Wc thus obtain 



E n = n(l - xi)q x + n{x\ - n)q + 2n(l - x 1 )^-x 1 q 1 + {x[ - n)q + 2n(x 1 - 



' {-Xiqi + (xi - n)q Q )(-nq ) 
-qi + (1 - n)q 

(26) 



This has the same ^/—nqo behavior as above. Setting 

qa = ~nr (27) 

we get up to order n 2 



E n = n(l - x 1 )q 1 - n 2 Xl r + 2n(l - x^-x^ - n 2 * 1 ^ + 2n 2 x 1 Jx~^ (28) 

y/-xiqir 

Optimization yields q\ = —x\ + 0(n), so 

E n = nx\(l — X\) — n 2 r + 2n 2 x ly /xir (29) 
Further optimization fixes x\ = 1/2, r = x\ = 1/8, which finally yields 

E n = -n + in 2 (30) 

The n 2 term is smaller that without replica symmetry breaking but still of wrong sign. We must allow for more 
breakings. 



4 



II. INFINITE ORDER REPLICA SYMMETRY BREAKING 

Now let us extend the above approach and consider k-step RSB. Let us consider the following representation for a 
matrix q: 

k 

q a f3 = (qd - 9fc)<W + - qi-ijE^ (31) 

8=0 

where q_i = 0, Xj < Xj_i and x k+1 = 1, Xq = n. It thus holds £( x i)£( x fc) = xj£( Xk ^ when Xj < x k . The eigenvalues 
of q a p are 

Afc = qd - qk 

k 

Aj = qd - Xi+iqi - ^ _ ^')* i = 0, • • • , fc - 1 

j=i+l 

A_i = nq + A (32) 
Translational invariance implies A_i = 0. Expressing q in the A's we get 

k 



i=0 

Its inverse has the same structure, so it can be written 



q a p = XkS a , + £ ^ (33) 



q a p = b a p = (b d - bk)S a p + - bi-i)£^ 

i=0 



1 yl/AM-l/Aj^, 
Xk ^ 

Let us define 

r Q/ 3 = {q~ r ) aa - {q~ 1 )ai3 (35) 

It has elements 

1 ^ 1 / 1 1 . 

A/. ^— '. £j Aj_i A, 

The replicated energy becomes 

A; k 

— = (xi - n)q Q + y(x j+ i - x i)q* + 2 — 1=^ + 2 7= — (37) 

n i=i v r o i=1 V r * 

where Xfe+i = 1. This can be rewritten as 

E n = nE™ + n 2 A n (38) 

with 

m = £{( — - -)Ai + 2(% '" Xi) } (39) 
i=l ^i+l v 7 ** 



with Xfe + i = 1 and 



A n = (l-i)^ + fc^ (40) 
xi n n n^/ro 



Only the terms in A„ depend explicitly on n. Using ro = T\ + (1/Ao — l/\i)/x\, which follows from (p6|), the 
dependence on Ao has been made explicit. The saddle point w.r.t. Ao yields 

A ° = n T , TTTJ « ri 2 xl(l - 3n 2 xl(x iri - ^)) + • •■ (41) 

(1 + Ao(xir! - ir )) J Ai 

in agreement with previous finding Ao = — nqo ~ +n 2 . We can now expand A„ in powers of n 

A n = - lUj - n 2 a;f (xiri - -!-) + n 3 xf (xin - -!-) H (42) 

Ai Ai 

With this form for A„, E n has to be optimized in x\, ■ ■ ■ , Xk and Ai • • • , A&. The n 2 term of E n has prefactor x\ , which 
is expected to vanish for fc — * oo. So infinite replica symmetry breaking cures the problem of negative variances. We 
have thus reached the main point of this work. It remains to be shown that the leading term (which is of order n 5 ), 
has a non-vanishing prefactor. Before doing that we calculate the average free energy. 

A. The average free energy 

Let us now calculate the 0(n) value of E n and later the leading correction, that will turn out to be of order n 5 . 
In the continuum limit one gets from eq. (B7p 



E (1) = / dxq{x) + 2 [ 
Jo Jo 



1 dx 

Tl^ (43) 
\/r(x) 



f 1 dy ... q'(x) 

\{x) J x y 2 X{y) X 2 (x) 

X(x) = -xq(x) + I dyq(y) — > A'(x) = -xq'(x) (45) 







We search a solution non-constant (7(2;) for < x < x, while 5(2;) = q for x < x < 1 for some a;. Variation with 
respect to q(z) yields 

1 f x dx 1 f z dx f x dy dx . . 

' ' ' =1 (46) 



r 3 / 2 (z)A 2 (z) J z xr 3 / 2 (x)X 2 (x) zX(z) 2 J Q r 3 / 2 (x) J z y 2 X 2 {y) ft r 3 / 2 (x) 

Differentiation with respect to 2 yields q'(z) = or 

1 2z f z dx 

2/ -T7TTT = (47) 



r 5 / 2 (z)A(z) r 3 / 2 (z) 7 r 3 / 2 (:r) 
Taking two more derivatives one obtains 



62 

Inserting this in eq. (fl7p one gets an equation for r alone. This yields 



r(z)X(z) = |- (48) 



x 6 6A 

with some A. Inserting this in eq. (0) one finds 



r(x) = 4 — » - ^ (49) 



G 

and finally from (flq) one obtains 



x = I (50) 



G 



= -^r- ( 51 ) 

K*) = (^) 

o3o9„5 

X(z) = — — (53) 



The average free energy is determined by 



= f + (1 - S)g + 2 f + 2^ = ^- (54) 

This value is approximately 4 times smaller than the exact value 7r/3. 

B. Cumulants of the free energy 

Next we consider the higher powers in n, coded in eq. (f|2|) for A„. For our purpose we may calculate it using the 
n = expressions for the Xi and A^. The terms of order rt°, n 1 and n 2 have too many powers of x\ and vanish for 
k — > oo, where xi — > 0. The n 3 term yields 

= x\r x ~ (55) 
Ai 

Using the continuum results ( |5l] ) it can be estimated by setting r\ ~ r{x\), Ai w A(xi) 



£ (5) - "^To (56) 



The correct approach is to calculate n and Ai at finite but large k. At n = this can be done quite simply. Taking 
derivatives of ( |39| ) with respect to Ai and x\ we obtain, 

1 1 z 2 -zi -1 _ Ai 2 _ 

a; 2 si r f/ J a; 2 Af xf yTi 



respectively. £2 can be eliminated, and we get 



l 

We donot need the value of X\ here, and obtain immediately 

„5 



A ' = ** r ' = IS? < 58) 



The prefactor in the estimate ([56]) has the same sign but is larger by a factor 1.32305. This result thus proves the 
leading behavior 



5! 7cr 
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(-P F )cum = ^L^- (60) 

which has the expected sign. In other words, in the variational approach free energy fluctuations are predicted to be 
of order 1/L 1 / 5 . 

III. SUMMARY 

In the present work we have extended the variational approach of Honeycutt and Thirumalai applied to directed 
polymers in 1+1 dimensional random media. This is a toy model for the more general field of interfaces in random 
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media, considered by Mezard and Parisi. In our problem a mapping on a Schrodinger equation is possible, which 
prevents the need of space-dependent propagators and the occurrence of wavenumbers. 

We have reconsidered the variational Ansatz of Honeycutt and Thirumalai. We have shown that it leads to a 
negative variance of the physical free energy, and thus is not a correct Ansatz. 

Next we have extended this Ansatz by allowing off-diagonal components of the matrix q a p. It was observed that 
an extra constraint is needed to make the problem well posed. The only natural and realistic candidate is spatial 
translational invariancc of the variational wavefunction in replica space. In doing so, we have seen that both the 
replica symmetric and one-step replica symmetry broken solutions still lead to a negative variance, though it has 
become smaller. We observe that more and more breakings improves the matter and that for infinite breaking the 
variance of the free energy has become zero (to leading order in the system size L). Infinite order breaking was also 
assumed by Mezard and Parisi. We find that the first non-vanishing cumulant is the fifth one, implying free energy 
fluctuations of order L x with \ — 1/5, which agrees with the estimate of Mezard and Parisi. Let us recall that this 
is the Flory or mean field estimate; loop effects should modify it to yield the exact result % = 1/3. 

This present situation of negative variances is reminiscent of the mean field model of spin glasses, SK-model, where 
an infinity of breaking was needed to obtain a vanishing zero-point entropy. Both effects are related with disappearence 
of unstable modes in the limit of infinite breaking. 
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